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1. Introduction 

The study of stochastic perturbations of dynamical systems is important in 
applications to realistic systems since their description contains uncertainties due to 
degrees of freedom that are not included, either because they are not measurable, or 
because their inclusion would complicate the analysis. A small stochastic perturbation 
can smooth the singular distributions inherent in many nonlinear dynamical systems. 
This makes it easier to study averages over such distributions as well defined 
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deterministic limits of smooth stochastic distributions. Besides being an ubiquitous 
fact of life - any dynamics in nature is stochastic to some degree, however weak - 
stochastic processes offer us a great freedom in picking problems and testing ideas. 

The properties of the weakly stochastic system are here determined from the 
unstable periodic orbits of the unperturbed (deterministic) system, decorated by the 
stochastic corrections. The central object in the theory, the trace of the evolution 
operator, is a discrete path integral, similar to those found in field theory and statistical 
mechanics. The weak noise perturbation theory, likewise, resembles perturbative field 
theory, and in the preceding paper [jlj] we developed such a perturbation theory for 
trace formulas for weakly stochastic chaotic dynamics in the standard field-theoretic 
language of Feynman diagrams. 

Here we approach the same problem from an althogether different direction; 
the key idea of flattening the neighborhood of a saddlepoint can be traced back to 
Poincare 0, and is perhaps not something that a field theorist would instinctively 
hark to as a method of computing perturbative corrections. In the Feynman 
diagram approach [0 we observed that the sums of diagrams simplify for saddlepoints 
corresponding to repeats of shorter periodic orbits, and were surprised by the 
compactness of the order er 2 correction. Here we explain this simplification in 
geometric terms that might be applicable to more general field theoretic problems. 

2. Stochastic evolution operator 

The periodic orbit theory allows us to calculate long time averages in a chaotic system 
as expansions in terms of the periodic orbits (cycles) of the system. The simplest 
example is provided by the Perron-Frobenius operator 

Cp(x') = j dx 5(f(x) — x')p(x) 

for a deterministic map f(x) which maps a density distribution p(x) forward in time. 
The periodic orbit theory relates the spectrum of this operator and its weighted 
evolution operator generalizations to the periodic orbits via trace formulas, dynam¬ 
ical zeta functions and spectral determinants & 0- For quantum mechanics the 
periodic orbit theory is exact on the semiclassical level ||, whereas the quintessentially 
quantum effects such as creeping, tunneling and diffraction have to be included 
as corrections. In particular, the higher order h corrections can be computed 
perturbatively by means of Feynman diagrammatic expansions (b) . Our purpose here 
is to develop the parallel theory for stochastic dynamics. In case at hand, already a 
discrete time 1-dimensional discrete Langevin equation |t], 0 , 

*£n +1 — f (*£rc) T 3 (1) 

with independent normalized random variables, suffices to reveal the structure of 
the perturbative corrections. 

We treat a chaotic system with weak external noise by replacing the deterministic 
evoluton ^-function kernel by the Fokker-Planck kernel corresponding to ( 0 ), a 
sharply peaked noise distribution function 

C a (x',x) = 6 a (f(x) - x'). ( 2 ) 

In ref. 10 we have treated the problem of computing the spectrum of this operator 
by standard field-theoretic Feynman diagram expansions. This time we formulate 
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the perturbative expansion in terms of smooth conjugacies and recursively evaluated 
derivatives. The procedure, which is relatively easy to automatize, enables us to go 
one order further in the perturbation theory, with much less computational effort than 
Feynman diagrammatic expansions would require. 

In the weak noise limit the kernel is sharply peaked, so it makes sense to expand 
it in terms of the Dirac delta function and its derivatives: 

_°° n rrm 2 3 

S*(y) = Y m I = S {y)+a 2 ^-S (2) {y)+a 3 —S (3) (y )+... .(3) 

z —' ml zb 

m —0 

where 

ak 

and the coefficients a m depend on the choice of the kernel. We have omitted the 
6^\y) term in the above because in our applications we shall impose the saddle-point 
condition, that is, we shift / by a constant to ensure that the noise peak corresponds 
to y = 0, so d CT (0) = 0. For example, if S a (y) is a Gaussian kernel, it can be expanded 
as 




co 

— cr zn 
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7T(J 
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s(y) 


n—0 
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-$ (2) {y) + ^<5 (4) (y) 


( 4 ) 


3. Stochastic trace formula 


We start our computation of the weak noise corrections to the spectrum of C a by 
calculating the trace of the nth iterate of the stochastic evolution operator C a for a one¬ 
dimensional analytic map fix) with additive noise cr. This trace is an n-dimensional 
integral on n points along a discrete periodic chain, so x becomes an n-vector x a with 
indices a, 6,... ranging from 0 to n—1 in a cyclic fashion 

/ n— 1 

[J dx a S a (y a ) 

a =0 

y a (x) = f(x a ) - Xa+l , Xn = X 0 . (5) 

To the order cr 3 the composition is simple: all compositions but one can be made, 
resulting in 

tr£” = J dx 5(f n (x) — x) + ~^ >j2 J dx\...dx n 

n — 1 

Y. ~ X l) ~ ■■^ 2 \f( x a) ~ x a+ 1) • • -8{f(xi) - X 2 ) + •• 

a —0 

n — 1 „ 

= tr£" + ^-cr 2 Y, / dx5 (2 \f n {x a ) - x a ) 

a—0 ' 
n — 1 p 

+ ^cr 3 Y / dx5 {3 \f n (x a ) - Xa) + 0(cr 4 ) . 


a—0 ' 


( 6 ) 
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At fourth order we get contributions from 8^\ as well as the two-point contributions 

n — 1 

‘>c = <"') + ^ 4 £ 

a —0 

+ 'T cr4 £ f dx a dx b 6 {2 \f J (xa) - x b )8 (2) (f k (x b ) - x a ) ,(7) 

a<b 

where j is the number of steps from points a to b on the cycle, and k is the number 
of step from b to a, so that j + k = n. 

If the map is smooth, the periodic points of given finite period n are isolated 
and the noise broadening cr sufficiently small so that they remain separated, the 
dominant contributions come from neighborhoods of periodic points; in the saddlepoint 
approximation the trace (|j) is given by 

trC—> £ (8) 

x c eFix/ n 

As traces are cyclic, e w ° is the same for all periodic points in a given cycle, independent 
of the choice of the starting point x c . Hence it is customary to rewrite this sum in 
terms of prime cycles and their repeats, 

OO 

tr K I saddles = £ n P £ ^ > ( 9 ) 

p r =1 

where p r labels the rth repeat of prime cycle p. 


J dx 8 {i) (f n (x a ) - x a ) 


4. Trace evaluated at a fixed point to all orders 


We now derive the perturbative expansion for a fixed point (n = 1) to all orders in 
cr. A fixed point and its repeats are of particular interest having the same interaction 
at every site, as does the usual field theory. What we do here is to formulate (and 
partially solve, in the sense of determing a few orders of the exact perturbation theory 
expansion) the field theory on finite periodic 1-dimensional discrete chains. 

Defining y = f(x) — x , we can write the fixed point trace as 

trC a = J dx 8 a (/ (x) — x) = J dy , * 8 a {y). (10) 

Expanding the kernel 8 a {y) as in (||) and integrating by parts, we see that all is well 
if we know the d/dx derivatives of l/y'{x). Replacing 

d dx d Id 
dy dy dx y'{x) dx 


we obtain in the saddlepoint approximation contributions to each fixed point of / 
evaluated recursively as derivatives of l/y'(x) 


/ 


dx 8^ (y) 


E (-!)* 

x:y(x )=0 


d k l 
dy k \y'{x)\ 



\y'(x)\ ' 
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The d/dy derivatives of 1 /y' are related to the d/dx derivatives of the map / by 
1 1 
,,r 


y 

d_l_ 

dy y' 
£_\_ 
dy 2 y' 
d 3 1 
dy 3 y' 


f'(.x) ~ 1 
f" 
(2/0 3 
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{y') 


'\7 


Here, and also later, we have relegated the fourth and fifth derivatives (A.l) to the 
Appendix. For example, for the second derivative of the delta function we have 


dxS^(y)= uTi 

x:y(x )=0 


/f \2 


y" 


, jy") 

’ (y'Y ( y'f 


where the sum is over all fixed points /(x) = x. In general, for n > 1 

d n 1 1 x ^ (2n — fci)! -[-j- f^ kl 

dy n \y'\ | y' 


(-1 r 


{fc^} 


{-y'Y 


i> 2 


(£\) k ‘k e \ ' 


( 12 ) 


(13) 


where is the £-th derivative and the sum is over all sets {kg} satisfying £ > 1, 
h > o, = n, and ^ £ki = 2n, that is, all partitions of 2n into n terms. The 
product contains only the finitely many nonzero kg for £ > 2. This formula may be 
proved in a staightforward manner by induction. It may be simplified substantially 
by defining a generating function 


T{z,x) 


5 > 


d ”- 1 i 
dy n ~ 1 y'(x) 


x:?/(a;)=0 


= exp 



f(n) z n d n 



(14) 


where = d/dy' acts only on y'. The sum in the exponential is formally a Taylor 
series expansion of 

T(z,x) = exp {— 2 _1 [f(x - d u ) + Ad u - x]} - , (15) 

u 

with u = y /z, f'(x) = A and /(x) = x. 


5. Smooth conjugacies 

The next step injects into Held theory a method standard in the analysis of fixed 
points and construction of normal forms for bifurcations, see refs. I]-®- The idea 
is to perform a smooth nonlinear change of coordinates that flattens out the vicinity 
of a fixed point and makes the map linear in an open neighborhood. This can be 
implemented only for an isolated nondegenerate fixed point (otherwise higher terms 
will contribute to the normal form expansion around the point), and only in a finite 
neigborhood of a point, as the conjugating function in general has a finite radius of 
convergence. Later we extend the method to periodic orbits, which are fixed points of 
the nth iterated map. 
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Let the fixed point of f{x) be x = 0 and the stability of that point be A = /'(0). If 
|A| 7 ^ 1, there exists a smooth conjugation h(x) satisfying h( 0) = 0 such that: 

f(x)=h{Ah~ 1 (x)). (16) 

In several dimensions, A is replaced by the Jacobian matrix, and one has to check that 
its eigenvalues are non-resonant, that is, there is no integer linear relation between 
their logarithms. If h(x) is a conjugation, so is any scaling h(ax) of the function for a 
real number a. Hence the value of h'( 0) is not determined by the functional equation; 
we shall set h'(0) = 1 . 

To compute the conjugation h we use the functional equation h(Ax) = f{h(x)) 
and the expansions 

f{x) =Ax + x 2 f 2 + x 3 f 3 + ... 

h{x) = x 4- x 2 h 2 + x 3 h 3 + ... . (17) 


In the present context absorbing the factorials into the definition of expansion 
coefficients turns out to be more convenient than the standard Taylor expansion. 
Equating recursively coefficients in expansions 

OO 

h(Au) - A h(u) = fm 

n =2 

OO OO / OO \ m 

^(A re - A)h n u n = Y, fmu m 1 + E ( 18 ) 


k—2 


yields 


h 2 = 


h 


h* = 


2/| + A(A - 1 )/ 3 


(19) 


2 A(A — 1) ’ 3 A 2 (A — 1)(A 2 — 1) ’ 

Noting that the left-hand-side of (|lj) generates denominators in ([ii]) are the same as 
those appearing in the Euler formula 


n (i+ tv ?) -1+j 


t 


t u 


3„.3 


+ 


Fu 


k=0 


— u (1 — w)(l — u 2 ) (1 — u)(l — rt 2 )(l — u 3 ) 




fc(fe-l) 

U 2 


k—0 


(1 — u) ■ ■ ■ (1 — u k ) 


we find it convenient to factorize h n as 
hn 


hn. — 


D n 

Computer algebra then yields 


Dr, = 1 - 


A 


1 - 


A 2 


\u\ < 1 , 


1- 


A" 


A 


(n + 2)(n-l) 


t >2 = /2 

63 = 2/ 2 2 +A(A-l)/ 3 (20) 

6 4 = (5 + A)/| - A(5 - 2A - 3A 2 )/ 2 / 3 + A 2 (A - 1)(A 2 - l)/ 4 . 


The expressions (A.2) for 65 and b§ are given in the Appendix. 
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Now that we have constructed the conjugation function for a fixed point, we turn to 
the problem of constructing it for periodic orbits. Each point around the cycle has 
a differently distorted neighborhood, with differing second and higher derivatives, so 
we need to compute a conjugation function h a at each cycle point x a ■ We expand the 
map / around each cycle point along the cycle, 

y a {4>) = fa(4>) - X a+1 = <pfa,l + 4?fa ,2 + ■ • ■ 
where x a is a point on the cycle, f a (4>) = f{x a + <j>) is centered on the deterministic 
orbit, and the index k in f a k refers to the fcth order in the (modified) Taylor expansion 

» „ 

For a periodic orbit the conjugation formula (|16|) generalizes to 

fM) = h a+1 {fMh- a i m 


at each point. The conjugationg functions h a are obtained in the same way as before, 
by equating coefficients of Taylor series, and assuming that the cycle stability is not 
marginal, |A| 7 ^ 1. The explicit expressions for h a in terms of / are obtained by 
iterating around the whole cycle, 

f n (x a +(/)) = ha^Ah,- 1 ((/)))+x a , ( 21 ) 


so each h a function is given by the derivatives given in the previous section acting 
on / n , evaluated at each cycle point a (recursive formulas for evaluation of such 
derivatives are given in Appendix A of ref. [lj). Again we have the freedom to set 
h' a (0) = 1 for all a. 

We shall also find it convenient to define partial stabilities along cycle laps 
f J {x a + 4>) = x b + h b {Ajh~ l (4>)) 
f k (x b + </>)= x a + h a (Akh^ 1 (</>)) 

where 

b—1 a —1 

Aj = f'(x c ), A k = f'(x c ) 

c—a c=b 

and AjA k = A. 

At this point we note that while in the deterministic case the cycle weight is given 
by the cycle stability, a quantity invariant under all smooth coordinate transfomations, 
the higher order corrections are not invariant. The noise is defined with respect to a 
particular coordinate system, and it has no invariant meaning. 


6. Trace formula for repeats of periodic orbits 

What is gained by rewriting the perturbation expansion for / in terms of the equally 
messy perturbation expansion for the conjugacy function hi Once the neighborhood 
of a fixed point is linearized, the repeats of it are trivialized; from the conjugation 
formula © one can compute the derivatives of a function composed with itself r 
times: 

f r (x) = h(A r h-\x)). 

One can already discern the form of the expansion for arbitrary repeats; the answer 
will depend on the conjugacy function h{x) computed for a single repeat, and all the 
dependence on the repeat number will be carried by factors polynomial in A r , a result 
that emerged as a surprise in the Feynman diagrammatic approach |[I|] . 
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The above observation enables us to move from a field theory constructed at a single 
point to a field theory which is translationally invariant on a periodic chain of arbitrary 
length r. The first and the second derivatives are 

4-f r (x ) = ti(A r h~ 1 (x))A r .. 

dx J w v v " h'{h.- l (x)) 

d 2 r( , A 2 r h"(A r h- 1 {x)) A r h"(h.-\x))h'(A r h-\x)) 

dx 2j [ ’ h'{h~ l {x )) 2 h'{h~ l {x)f ' 1 ’ 

The third derivative is too long to write down here. Evaluation at the fixed point 
x = 0 yields 


—ffO) = A r (A r -l)h^ 
f r ( 0) = A r (A 2r - 1 )h {3) - 3A r (A r - 1 )h {2)2 ■■■ . 


-f/ r ( 0) =A'\ 

ax 
d?_ 
dx 3 


(23) 


(For brevity we shall often replace M r )(0) —> = r\h r , /^(O) —> f^ = r\f r 

throughout, where the superscript ^ labels the rth coefficient in the standard Taylor 
expansion.) In general, the fcth derivative of a function composed with itself r times 
is a polynomial in A r of degree k 

it- fc 


dx k 


/ r (0) = 5>A*\ 


.7=1 


The coefficients Cj are in turn expressed in terms of the first k derivatives h^> of the 
conjugation function. 

Now we can reevaluate the fixed-point derivatives of sect. |] for the case of rth 
repeat of a fixed point 

y(x) = f npr {x) - x 

in terms of the conjugating function h by substituting (^) into &■ 


(24) 


1 d 2 


3! dy 2 y'(x) 


A r (1 + A r ) / , 2 x 
(A r - l ) 3 ^ 2 


1 d 3 


1 


4! 9y 3 y'(x) 


c A r (A r + l ) 2 3 r 5A 2r + 8 A r + 5 
= - 5A —- TTjh 2 + A - TTZ - 7 T 7 i-«2ft3 


(A r - l ) 4 


(A r - l ) 4 
r A 2r + A 1 " + 1 

— A -—-—- tl4 

(A r - l ) 4 

The 4th and 5th order derivatives (A.3) are relegated to the Appendix. 


(25) 


(26) 


6.2. Repeats of periodic orbits 

The above fixed-point formulas carry over to the case of general periodic orbits, if / 
is identified with the ?rth iterated map at a particular point along the cycle a from 
which we define h a as in (Ell) . The second derivative required for the a 2 correction 
follows immediately, with ( p5| ) evaluated for each cycle point in the prime cycle trace 
(§• It is easily checked that this second derivative (|2^) has exactly the same form 
as what we obtained in a much more involved manner in ref. [0 . The somewhat 
more complicated sum of Feynman diagrams is here replaced by (25), the conjugation 
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function determined from the iterated map, with h carrying all dependence on the 
higher derivatives of the original map. 

The factorisation is not quite as simple for higher orders in the trace formula. 
For a general non-Gaussian case at order er 3 we have three terms from (p 6 |), each 
of which may be resummed separately. The calculation is entirely analogous to the 
second order calculation in ref. |]I| . Powers of A r are moved around in order to express 
( p 6 | ) in terms of |A| r < 1 and generate convergent geometric series; each of the three 
terms is expanded using identities such as 


(1 + *) 2 
(1 - x ) 4 


Y *(1 + 2fc 2 )x fc , 

k =0 


(27) 


so that repeats r can be resummed. (The remaining identities ( |A. 6 | )-( [A. 8 | ) used here 
are relegated to the Appendix.) This yields the saddlepoint approximation to the 
spectral determinant in product form up to order er 3 


OO 

det(l - zC a ) = J|(l ^ t Pt k) 

p k =0 


(28) 


fp,* = Py XF ex P (\ a *> 2w pl + pi“3+ 0(<r‘)j 

■S=(HifEK f -‘ > 3 >) 

a 

y; + a +1 )( 5+si + y 

a a 

(1 + k )(2 + 2k + fc 2 ) > ^( 4 -) 

2 

a 

The meaning of this factorization is that in the saddlepoint approximation the spectral 
determinant is composed of local spectra evaluated for each prime cycle separately, 
with the index k labelling the £;th local eigenvalue. 

From er 4 onward further terms come into play, as we now describe. 


7. Two-point integrals 


In general, a saddlepoint at each cycle point may be expanded in second and higher 
derivatives (|j), with the trace to nth order receiving a single <5^ integral contribution 
as in (|Gj). At fourth order we need to also include the two-point integral contribution 
(yi). Setting yj = fi{x a ) — Xb , yu = f k (xb) — x a , and integrating by parts, we obtain 
for each pair of points along the cycle 


dyjdy k 


| f j '(x a )f k '{xb) - 1 


■S^(yj)6 (2) {y k ) 


d 2 d 2 1 

dyj dy 2 k \f 3 '(x a )f k '{x b ) - 1 | ' 


This expression looks innocent enough, but its evaluation requires attention to some 
subtle points, and more algebra than a human would be willing to handle. Performing 
the required derivatives taxes the ability of general purpose symbolic algebra packages, 
so the results in this section are obtained with a program written in C. 
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Consider evaluating partial derivatives of a function q with respect to, say y :j 
keeping y constant in terms of x a and x b derivatives. The solution is to consider an 
arbitrary infinitesimal transformation which respects the dyk = 0 constraint: 


dq =P^dx. 


dyj = —2-dx 


dyk = dx. 


dx a 

dyj 

8 X a 

dyk 

dXn 


dq 

a— dx b 
dx b 

dyj, , 
w —dxb 
dx b 

dyk , n 

-— dx b = 0 

OXb 


dg dy k 
dx a dx b 


dq dyk 

dx b dx a 


dyj 

8Xn 


dyk 

dx b 


dyj 

dx b 


dyk 
dx a 


These equations can then be solved for the ratio of dq and dyj , 
dq 
dyj 

and similarly for dq/dyk- 

Sixth order corrections require a three point integral, which we leave as an exercise 
for the reader. 

If for some reason the noise kernel contains a linear cqcr term in the expansion (^]), 
for example if the evolution operator is weighted, and the shift that this causes to the 
saddlepoint is not taken into account, the relevant second order correction includes 
the term 

_d _ d _1_ 

dyj dyk \(dfj(x a )/dx a ){df k (x b )/dx b ) - 1| 

= [( A + 1)^ 2) >4 2) + 2A ^i 2)2 + 2A ^i 2)2 - A ^i 3) - 

This we have written down just to illustrate the form of a 2-point term. In our 
application the saddlepoint condit ion sets oi = 0, and what is really required is the 
two-point fourth order corr ectio n (A.4) and (in case of non-Gaussian noise) the two- 
point fifth order correction ( A.5 ). The corresponding expressions are too cumbersome 
for the main text, so we relegate them to Appendix. 

There are som e co nvenient cross-checks on the algebra. The results reduce to 
those of @) and ( A.3 ) for single repeat r = 1, by setting h a = h b = h, A ? = A, 
A*, = 1, and remebering that in order to retain integer coefficients, we have here 
replaced power series coefficients h r by the Taylor expansion derivatives h ^ = r\h r . 
This reduction to the one-point case also means that for the Gaussian case (02 = 1, 
04 = 3), expression (Jrj) can be written as an unrestricted double sum, 

^ v- ^ d 2 d 2 1 

R 2-^ 


E 

x(=Fix/ n ab 


dy] dy\ \(df^{x a )/dx a )(df k (x b )/dxb) - 1| 


where it is understood that j = n and k = 0 when a = b. 

Now the stage is set to return to the investigations of ref. jl] , where the a 2 term 
was evaluated anlytically, while the <r 4 term was only estimated from the numerically 
computed leading eigenvalue. 


8. Numerical tests 

Here we continue the calculations of Sect. 5 of ref. Jl|], where more details and 
discussion may be found. We test our perturbative expansion on the repeller of the 1- 
dimensional map f(x) = 20[(1 /2) 4 — ((1/2) — x) 4 ]. This repeller is a nice example of an 
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“Axiom A” expanding system of bounded nonlinearity and complete binary symbolic 
dynamics, for which the deterministic evolution operator eigenvalues converge super- 
exponentially with the cycle length. We compute the leading eigenvalue of the 
evolution operator (the repeller escape rate) in the presence of Gaussian noise, using 
two complementary approaches. The perturbative result in terms of periodic orbits 
and the weak noise corrections is compared to a numerical eigenvalue obtained in 
ref. jlj] by approximating the dynamics by a finite matrix. In the preceding paper we 
compared the numerical eigenvalue with the a 2 result and estimated the coefficient of 
cr 4 to be approximately 38. Here we compute the order cr 4 coefficient to 14 digits 
accuracy, and estimate the cr 6 term. In the Feynman diagram language, the cr 4 
contribution is a “2-loop” calculation, albeit one of relatively simple kind where space 
integrals are replaced by discrete sums over a finite periodic chains. The conjugation 
functions ©) at each point around the cycle are obtained by a recur sive evaluation 
of (|l^), and then substituted into ( |25| ) for the second order and (|Aj) for the fourth 
order. The expansion for the spectral determinant is obtained by differentiating 
lndet(l — zC a ) = trln(l — zC a ) with respect to z, multiplying through by the 
determinant, and equating coefficients order by order in 2 and cr, as in ref. (l). That 
is, we define coefficients 


tr Of = Y, ( 'r,.0 J , det(l - zC a ) = 1-^12 

j—0 n=l j=0 

and obtain the cumulants Q n ,j recursively as 
1 


,J Z ° > 


Qn, 0 — 

n 


C n , 0 — E Qk, oC n —k,Q 


k= 1 


Qn, 2 — 

n 


n— 1 


C n , 2 - E ( Qk^Cn—kfi T Qk,oC'n—k,2') 


k=l 


Qn,A — 


n— 1 


Cn,4 — 'y ] {Qk,A^ri-k,0 + Qk,2C n —k,2 + QkfiCn—k,A) 


k=l 


where it is understood that the sums do not contribute when n = 1. 

The noiseless, zeroth order eigenvalue equation det(l — zCq) = 0 is solved by 
Newton’s method to find the leading eigenvalue vq = z~ 1 at a = 0, and the higher 
order equations give the noise corrections v{cr) = i'o + v 2 (j 2 + iqcr 4 + 0(cr 6 ) in terms of 
z/o and the expansion of the determinant. Expanding the spectral determinant order 
by order in 2 and cr we find 

det(l - zLa) = F— F w z - F 02 (t 2 - F 20 z 2 - Fi 2 za 2 - F 04 cr 4 - .. .(29) 


where the expansion coefficients are 


77 > 1 \ A Qm ,0 -r—\ \ ' ^Qm,0 \ ^ Qm ,2 

^ ~ 4 / y > ** 10 — / y m — 1 > ^ 02 — / y ,.m 

m=l m—1 m=l ° 


f 20 = 


m(m - l)Q m ,o 


p _ SO m Qm ,2 p _ V 

2v™~ 2 ’ 12 E 2 ’ F ° 4 E 

u m—1 u 


Qr 


=1 


are obtained from derivatives of (^). Finally we expand v = 2 1 in powers of cr 2 and 
equate coefficients of powers of cr to obtain 


v 2 = 


Fo 2 v 0 _ F 2 qFq 2 — 2F12F10F02 + F 04 F 10 + F 10 F 02 u 0 2 


F10 


v 4 = 
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0.02 0.03 0.05 0.07 0.1 



Figure 1. Numerically computed eigenvalue minus known terms, that is, 
v(cr) — vo — cr 2 V 2 — cr^v^ (dots), together with estimated next term 2700cr 6 (solid 
line). Deviations occur at small a due to discretization errors in the numerical 
eigenvalue of a few times 10“ 7 , and at large cr due to further omitted terms 

(o' 8 ---)- 


n 

*4) 

Vi 

V ‘1 

i 

0.308 

0.42 

2.2 

2 

0.37140 

1.422 

32.97 

3 

0.3711096 

1.43555 

36.326 

4 

0.371110995255 

1.435811262 

36.3583777 

5 

0.371110995234863 

1.43581124819737 

36.35837123374 

6 

0.371110995234863 

1.43581124819749 

36.358371233836 


Table 1. Significant digits of the leading deterministic eigenvalue and its a 2 and 
cr 4 coefficients, calculated from the spectral determinant as a function of the cycle 
truncation length n. Note the super-exponential convergence of all coefficients 
(the n = 6 result is here limited by the machine precision). 


The perturbative corrections to the leading eigenvalue (escape rate) of the weak- 
noise evolution operator are given in table [j], showing super-exponential convergence 
with the truncation cycle length n. The super-exponential convergence has been 
proven for the deterministic, v 0 part of the eigenvalue |l9|, |o), but has not been 
studied for noisy kernels. It is seen that a good first approximation is obtained already 
at n = 2, using only 3 prime cycles, and n = 6 (23 prime cycles in all) is in this 
example sufficient to exhaust the limits of double precision arithmetic. The exact 
value of M 4 = 36.358... is encouragingly close to our previous numerical estimate [Jlj] 
of 38. As in the preceding paper jlj, we subtract the known terms in the expansion 
from the numerically evaluated eigenvalue, and obtain a good fit to the next term, 
approximately 2700cr 6 , see fig. [l[ 






Stochastic trace formulas 

9. Summary and outlook 


13 


We have formulated a perturbation theory of stochastic trace formulas based on 
smooth nonlinear field transformations around infinitely many chaotic saddle points 
(unstable periodic orbits). In contrast to previous perturbative expansions around 
vacua and instanton solutions, the location and local properties of each saddlepoint 
must be found numerically. In addition, every interaction term depends on the 
position, hence also the classical periodic solution at which it is evaluated. 

Even though in the model calculation we have chosen a Gaussian one-step noise, 
the accumulated noise along a trajectory is distorted by the nonlinear flow, and orbit- 
by-orbit the noise corrections are not Gaussian. Gaussian noise has thus no privilaged 
role in nonlinear dynamical systems. 

The key idea is this: Instead of separating the action into a quadratic and 
“interaction” parts, we first perform a nonlinear held transformation (“smooth 
conjugation”) which turns the saddle point into an exact quadratic form. The price 
one pays for this is the Jacobian of the nonlinear held transformation - but it turns out 
that the perturbation expansion of this Jacobian in terms of the conjugating function is 
order-by-order considerably more compact than the Feyman-diagrammatic expansion. 

We have resummed repeats of prime cycles to third order in the noise strength, 
carried out the numerical tests to fourth order, given the trace formula for general 
periodic orbits to hfth order, and for a hxed point to all orders. 

The rapid rate of increase of the numerical coefficients conhrms the expectation 
that the series is asymptotic, and is to be used with caution, as long as no sophisticated 
summation beyond all orders is implemented. 

The smooth conjugacy method of perturbation expansions can be extended to the 
case of higher dimensions and continuous time dynamics (stochastic flows), but the 
main interest comes from the observation that we have a new method of evaluating 
perturbative corrections to saddlepoints of path integrals. In quantum mechanics 
and field theory the perturbative corrections do matter, and the method might have 
applications there, in particular to the Ti expansion of semiclassical periodic orbit 
theory. If efficient methods are found for computing numerical periodic solutions of 
spatially extended systems, the method might apply to the field theory as well. 


Appendix A. Appendix: Some algebra 


We collect here some of the formulas used in our calculations, but too long for the 
patience of a casual reader. 

Evaluation of the trace of a fixed point ( 0 ) continued to 4tli and 5th order: 


d A 

1 

/ (5) 

15/"/"" 

10(/'") 2 

105(/") 2 /"' 

105(/") 4 

dy A 

y' 

{y'Y 

{y'Y 

{y'Y 

{y'Y 

{y'Y 

d 5 

1 

/ (6) , 

2i/"/(5) 

35/,/////// 

210 {f") 2 f"" 


dy 5 

y' 

0 v'Y 

{y'Y 

+ {y'Y 

{y'Y 




280 f"{f") 2 , 1260(/") 3 /'" 
{y'Y + {y'Y 0 

945 (/") 5 

{y'Y 1 ' 

(A.l) 


Perturbative formulas for 5th and 6th order conjugacy coefficients, continuation 
of @): 


fo 5 = 2(7 + 3A + 2A 2 )/| + A (-21 + A + 6A 2 + 11A 3 + 3A 4 ) / 2 2 / 3 
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+ 2A 2 (A - l) 2 (3 + 5A + 4A 2 + 2A 3 ) / 2 / 4 + 3A 2 (A - 1)(A 3 - 1 )/ 3 
+ A 3 (A — 1)(A 2 — 1)(A 3 — l)/ 5 
b 6 = 2 (21 + 14A + 14A 2 + 8A 3 + 3A 4 ) /f 

+ A (-84 - 21A - 6A 2 + 25A 3 + 44A 4 + 27A 5 + 14A 6 + A 7 ) /f/ 3 
+ 2A 2 (14 + A - 7A 2 - 9A 3 - 11A 4 + A 5 + A 6 + 7A 7 + 3A 8 ) /|/ 4 
+ (-1 + A) 2 A 2 / 2 [(28 + 43A + 50A 2 + 43A 3 + 22A 4 + 6A 5 ) /f 
+A (-7 - 12A - 10A 2 - 3A 3 + 7A 4 + 10A 5 + 10A 6 + 5A 7 ) f 5 ] 

+ (-1 + A) 3 A 3 (l + A + A 2 + A 3 ) 

x [(7 + 7A + 4A 2 ) / 3 / 4 + A (-1 - A + A 3 + A 4 ) / 6 ] . (A.2) 

Continuation of the derivative evaluation (|26|) - the 4th and 5th order derivatives: 


1 <9 4 1 

5! dy 4 y'{x) 


<V3) 

- 3 (7 + 10A r + 7A 2r ) h\h 3 + 3 (l + A r + A 2r ) h 2 3 
+2 (3 + 2A r + 3A 2r ) h 2 h A - (l + A 2r ) h s } 


_1_ 

6! dy 5 y'(:r) 


- 28(1 + A r ) 2 (3 + 4A r + 3A 2r ) hlh 3 
+ 14(1 + A r ) 2 (2 + A r + 2A 2r ) h 2 h 4 

- (7 + 14A r + 18A 2r + 14A 3r + 7A 4r ) h 3 h 4 
+ 7 (4 + HA r + 15A 2r + HA 3r + 4A 4r ) h 2 h 2 3 

- (j + 12A r + 12A 2r + 12A 3r + 7A 4r ) h 2 h 5 
+ (l + A r + A 2r + A 3r + A 4r )h 6 } . 


The two-point fourth order correction (continuation of calculations of sect. 


0: 


d 2 d 2 


1 


dy 2 dy 2 \(dfi{x a )/dx a )(df k (x b )/dx b ) - 1| 
A 


cycle 


(A-l) 5 


{(1 + A) [A 2 fc (60- 72+ Uh^h^ + 9h^ 2 - 


+A 2 (60/if )4 - 72+ 12/4 2) /4 4) + 9 
(3 + 10A + 3A 2 ) [A k h ( b 2) (l2h^ 3 - 9 + hty 
+A jhP (l2 h™ 3 - 9 h[ 2) h[ 3) + h^)] 

(1 + A)(l + 10A + A 2 ) (3 h {2)2 - ( 3 h 


(A.4) 


( 2 V _ . ( 3 ) 
b n b 


)}■ 


and the two-point fifth order correction is: 


d 3 d 2 1 

dy 3 dvl \(dp(x a )/dxa)(df k (x b )/dx b ) - 1| cycle 

(A ^ 1)6 {(1 + A)A^ (-360fti 2)5 + 600/7.£ 2)3 /i£ 3) - 180h^h^ 2 
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-120fti 2)2 fti 4) + 30/i,i 3) /ii 4) + 15/ii 2) /ii 5) - , 

+ 6(1 + 3A + A 2 )A 2 k h[ 2) (-60fei 2)4 + 72h (2) \ (3) - 9 h ( F - 12fc< a) fc< 4) + h i 5) ) 


+ 


15(1 + 7A + 7A 2 + A 3 )h (2Y - 2(2 + 13A + 13A 2 + 2A 3 )h (3) 


Afe (-mi 2)3 +9h^h {3) -/i< 4) ) 

15(1 + 18A + 42A 2 + 18 A 3 + A 4 5 )/if )3 


(A.5) 


-2(5 + 82 A + 186 A 2 + 82A 3 + 5A 4 )h (2) h (3) 

+(1 + 14A + 30A 2 + 14A 3 + A 4 )h ( b 4) (-3 h ( F + h (3) ^j 


+ 3 


-30(3 + 17A + 17A 2 + 3A 3 )h (2) + 5(19 + 101A + 101 A 2 + 19A' i )/i^ 2r h ( b 
-10(1 + 5A + 5A 2 + A 3 )h (3)2 

-2(7 + 33A + 33A 2 + 7A 3 )h (2) h { f + (1 + 4A + 4A 2 + A 3 )^ 6 7 8 
+ -630(1 + 2A + A 2 )hf )b + 30(31 + 58A + 31A 2 )h (2)3 h (3) 

-60(4 + 7A + 4A 2 )h (2) h (3)2 - 15(11 + 18A + 11A 2 )h (2)2 h { b 4) 

+2(17 + 26A + 17A 2 )h (3) h (4) + 6(3 + 4A + 3A 2 )/if - (1 + A + A 2 )h { b ] A 2 ) . 

The repetition number dependent prefactors in ( |2fi| ) are turned into power series in 
A rk using identities 


^,,( 2 )^( 3 ) 


A jh { a ] 


1 + x 

(i -xY 


= E( fc+1 ) 


2 x k 


fc =0 

5 + 8a: + 5x 2 


(1-x) 4 

1 + X + X 2 

(1-z) 4 


E (1 + *) (5 + 6 k + 3 k 2 ) x k 


k> 0 


— t; E 0 + k) (2 + 2 k + k 2 ) 


(A.6) 

(A.7) 

(A.8) 


fc >0 
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